Networking Science manuscript No. 

(will be inserted by the editor) 



Path Selection for Quantum Repeater Networks 

Rodney Van Meter • Takahiko Satoh • Thaddeus D. Ladd • 
William J. Munro • Kae Nemoto 



(N 



O 

in 

(N 



Or 
i 

G 



> 

in 
in 

in 

O 

(N 



Received: date / Accepted: date 

Abstract Quantum networks will support long-distance quan- 
tum key distribution (QKD) and distributed quantum com- 
putation, and are an active area of both experimental and 
theoretical research. Here, we present an analysis of topo- 
logically complex networks of quantum repeaters composed 
of heterogeneous links. Quantum networks have fundamen- 
tal behavioral differences from classical networks; the deli- 
cacy of quantum states makes a practical path selection al- 
gorithm imperative, but classical notions of resource utiliza- 
tion are not directly applicable, rendering known path se- 
lection mechanisms inadequate. To adapt Dijkstra's algo- 
rithm for quantum repeater networks that generate entan- 
gled Bell pairs, we quantify the key differences and define 
a link cost metric, seconds per Bell pair of a particular fi- 
delity, where a single Bell pair is the resource consumed 
to perform one quantum teleportation. Simulations that in- 
clude both the physical interactions and the extensive classi- 
cal messaging confirm that Dijkstra's algorithm works well 
in a quantum context. Simulating about three hundred het- 
erogeneous paths, comparing our path cost and the total work 
along the path gives a coefficient of determination of 0.88 or 
better. 
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1 Introduction 

A routing algorithm chooses a path on a graph, and consists 
of two parts: a definition for the cost of a single link, and a 
function for calculating the cost of a path based on those link 
costs, allowing us to extend a single point-to-point channel 
to a richer network. Dijkstra's Shortest Path First algorithm, 
for example, takes a simple scalar cost for each link and 
treats the sum of link costs as a cost for a candidate path JT|. 
The emerging field of quantum communication has, to date, 
experimentally demonstrated the basic principles of entan- 
gled quantum networking |2 3 4 5 6|. and laid the theoret- 
ical foundations of creating long-distance, high-quality en- 
tanglement J7]|8]|9l, but topologically has considered primar- 
ily channels and linear networks, leaving us with an urgent 
need for a path selection mechanism as quantum networks 
develop. 

Quantum key distribution (QKD) is probably the most 
prominent use of quantum communication, and commercial 
products are available II10II1 llfPH . In QKD, quantum effects 
(and a large dose of classical statistics) are used to detect the 
presence of an eavesdropper on the channel. QKD generates 
streams of shared, secret, random bits, which can then be 
used to key a cryptographic session, such as an IPsec tun- 
nel lll3lll4|[T5l[T6l . To extend QKD over distances exceed- 
ing a certain limit (a few hundred km in telecom fiber today) 
or to create networks without direct links between all pairs 
of nodes, we can use trusted relay nodes and direct optical 
switching. Such metropolitan-area QKD networks have al- 
ready been demonstrated in Boston, Vienna, Geneva, Tokyo, 
Hefei, and other cities lflT1[T3l[T6l[T7l . 

Long-distance QKD without the trusted relay nodes, as 
well as most other distributed quantum applications, can be 
achieved through the use of entanglement MT8TT9I (entan- 
glement and other italicized vocabulary will be explained 
in more depth in Sec. |2j. Physical links can create entan- 
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glement over short distances; in multi-hop networks, many 
short-distance entangled states are used to forge the end-to- 
end entangled states consumed by applications. Complex, ir- 
regular network topologies with hops of varying length and 
quality, from quantum LANs through MANs to WANs, will 
result as quantum systems grow. The devices that will bridge 
long distances and route connections through networks are 
called quantum repeaters {8). 

Entangled states have uses beyond QKD. The quantum 
communications toolkit includes such building blocks as quan- 
tum Byzantine agreement, distributed leader election, secure 
multi-party communication, and distributed arithmetic for 
building large-scale quantum applications ll20ll2Tll22ll23l . 

Thus, heterogeneous quantum networks hover in our fu- 
ture, but substantial theoretical and engineering hurdles re- 
main to be cleared. The first, theoretical studies on entangled 
quantum networks focused primarily on an abstract model 
consisting of a linear chain of repeaters, with a power of two 
number of hops of identical length and quality [8,24|. Re- 
cent work [25 26 27 1 has targeted more realistic chains of 
repeaters, relaxing those constraints. Here, we analyze the 
behavior of more complex network topologies, as in Fig. Q] 
In a network of heterogeneous links and irregular topology, 
path selection affects both the performance of individual 
connections and global network load. Purify-and-swap quan- 
tum repeaters are a theoretically well-developed approach, 
and are strong candidates for field deployment, and we use 
them as our subject. 




a route: A =5 D => E => C AAO 

P route: A => B => E => C xOO 

v route: A => B => C x x 
V J 



Fig. 1 Path selection is a critical problem in all networks; in quantum 
networks, due to the delicacy of quantum information, it may deter- 
mine whether or not a connection can be made successfully. 



In this paper, we apply Dijkstra's algorithm for ranking 
candidate paths to quantum repeater networks JT]. We eval- 
uate our algorithm for its ability to select a path that both 
maximizes the throughput of end-to-end connections and 
minimizes global work for this class of repeater. Detailed 
simulations of both the physical interactions and the clas- 



sical messaging confirm reasonable agreement between the 
calculated path cost and the expected throughput. 

Our proposed link cost is the inverse of the throughput 
of the link, measured in Bell pairs per second of a particu- 
lar fidelity. Other candidates for link cost, including lower- 
level metrics such as the number of laser pulses and quantum 
measurement operations, are found to be useful for evaluat- 
ing the total work actually consumed on a path, but are poor 
metrics for prioritizing a link because they reflect the phys- 
ical link characteristics but not the system factors that are 
equally important influences on end-to-end performance. 

We present the results of two sets of simulations of var- 
ious paths using four different qualities of links. The first 
set of forty-six paths vary in length from one to nine hops, 
while the second set covers 256 link combinations in four- 
hop paths. Across both data sets, the coefficient of determi- 
nation is 0.88 or better between the path cost and the total 
work performed (counted as the number of quantum mea- 
surements performed along the whole path), supporting our 
choice of link cost and the effectiveness of Dijkstra for this 
type of quantum network. Comparing the results of pairs 
of simulations, the path with the lower cost also has higher 
throughput in more than 80% of all tested cases. We demon- 
strate that, in direct analogy to classical networks, the per- 
formance of a quantum path will be limited by the through- 
put of the bottleneck link, while total work is a function of 
both the path length and the quality of all the links. 

To build a complete argument, we discuss the differ- 
ences between quantum and classical networks and the diffi- 
culties encountered (Sec. [2j. After defining the problems in 
path selection and proposing several solutions (Sec. [3), we 
evaluate those solutions via simulations that answer a series 
of specific questions about the behavior of quantum repeater 
networks (Sec.|4). 

2 Difficulties: Differences Between Quantum and 
Classical Networks 

So far in this paper, we have discussed quantum commu- 
nication technologies and outlined where we apply classi- 
cal techniques to the quantum problems. However, we have 
not specifically articulated the fundamental differences that 
make the merger of classical and quantum networking con- 
cepts less than straightforward. There are both theoretical 
results and practical reasons for believing that the answers 
to our questions require thought, rather than simply assert- 
ing that classical and quantum networks can use the same 
solutions. 

The differences stem from several sources: the engineer- 
ing difficulties of creating and protecting quantum states and 
the real-time decay of quantum information; the impact of 
(probabilistic) photon loss; and the fundamentally proba- 
bilistic nature of some quantum operations. These issues 
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manifest themselves both locally and globally, requiring in 
some cases additional classical messages to be exchanged, 
further exacerbating the problems. 

Perhaps most importantly, we cannot make independent 
copies of quantum information, a fact known as the no cloning 
theorem J28[. In classical systems, data is copied as it is 
transmitted; if the data is lost, retransmission is done with 
no harm beyond a performance penalty. With only a single 
copy of quantum data, we are forced to treat that data with 
extra care. Generally, quantum networks build generic states 
which are then used to teleport precious data Q. To under- 
stand teleportation, we must first discuss entanglement; after 
teleportation and its uses, we can complete the role of quan- 
tum repeater by discussing imperfect quantum states. 

2. 1 Entangled States 

Quantum communication involves the sending and receiving 
of quantum information, encoded in either the state of an in- 
dividual photon, or a quantum state of a stronger optical or 
microwave pulse. The state may be sent through a waveg- 
uide such as an optical fiber, or through free space, even via 
satellite l5l fT2l[T6ll29ll30ll3T1 . A transmitted state may stand 
alone, or be part of a larger, entangled quantum state. 

Entanglement, which Einstein referred to as "spooky ac- 
tion at a distance", is the condition of a quantum state in 
which the states of subsystems are not independent: opera- 
tions on one part of the system can affect the state of other 
parts, regardless of physical distance. This effect does not vi- 
olate the principle of relativity, because it cannot be used to 
transmit information instantaneously; interpreting the mean- 
ing of the distant parts requires the use of classical informa- 
tion acquired during the local quantum operations. 

That classical information about the quantum states is 
acquired via measurement of the states. Measuring a single 
quantum bit (qubit) produces one classical bit whose value 
depends on the details of the quantum state, and destroys 
any entanglement of the qubit with others. The actual mea- 
surement mechanism is of course technology-dependent, but 
often involves interferometric optical setups with very sensi- 
tive photodetectors, allowing polarization or other character- 
istics of an optical state to be determined, detecting a single 
photon or measuring the direction of spin of a single elec- 
tron. For the quantum states of interest here, the classical 
information gained is guaranteed by quantum mechanics to 
be random, and must be transmitted to the remote site, which 
of course can only be done at the speed of light. 

The particular type of entanglement we require is called 
a Bell pair. Although entanglement is location independent, 
when we refer to a Bell pair in this paper, we mean a ge- 
ographically distributed Bell pair consisting of one qubit at 
one location entangled with a qubit at another location. En- 
tanglement comes in many forms and can involve more than 



two parties, but in this paper we focus on the two-party Bell 
pairs both as intermediate resources and the end goal of the 
communication session delivered to the application layer. 

Quantum communication systems use entanglement in 
various ways, and for different purposes. Some uses require 
that the end nodes store quantum data in memory for ex- 
tended periods of time (many round-trip times); others uti- 
lize quantum states of light, but measure the state directly 
upon receipt without storing, converting the quantum state 
to a classical value. One of the greatest challenges in quan- 
tum communication is dealing with memory lifetimes that 
are low compared to round-trip latencies. 

2.2 Teleportation and Entanglement Swapping 

To teleport a qubit, we begin with a Bell pair with one end 
held at our source and the other held at the destination. Lo- 
cal entangling operations are performed on the data qubit 
and the source-end Bell pair qubit, then both are measured, 
giving two classical bits of information. Those two classical 
bits are transmitted to the destination, which then performs 
local quantum operations on its Bell-pair qubit, dependent 
on the values of those bits. This recreates the original data 
qubit at the destination. 

Entanglement swapping (Fig.|2j splices two short-distance 
Bell pairs into one longer-distance Bell pair: in Fig. [2] (A ■<-> 
B) + (B o C) => (A -n- C) combines the two left gold 
arcs to form the green one above. Swapping is a form of 
teleportation; it can be viewed as using the (B •<-» C) Bell 
pair to teleport the B end of the (A <H> B) pair to C. Clas- 
sical information must be sent to both ends of the new en- 
tangled connection. Managed carefully, each swapping step 
can double the span of our entangled Bell pair. 




1st 2nd 1st 

Fig. 2 Symmetric nested entanglement swapping, splicing shorter- 
distance entangled pairs to create longer-distance entanglement, is ac- 
complished easily in n steps for 2™ hops of the same quality. The 
legend and arrows show the order in which entanglement swapping is 
performed. The arcs represent entanglement over different distances, 
created via the swapping of shorter entangled pairs. 



2.3 Imperfect Quantum States 

In real systems, Bell pairs are imperfect, and have a char- 
acteristic known as fidelity, which represents the accuracy 



4 



of our knowledge about the quantum state, with F = 1.0 
representing perfection. Mathematically, the fidelity is the 
amount of overlap between the state we intended to create 
and an average over an ensemble of states we actually cre- 
ated, either experimentally or in simulation, and is defined 
as F = (ty + \p\\l/ + ), where p is the density matrix represen- 
tation of the state we have, and is the Bell state we 
are trying to create. In experiments, fidelity is determined 
statistically using large numbers of repetitions, creating a 
quantum state and measuring it to see if the state was prop- 
erly made and stored. In simulations, various imperfections 
are numerically introduced on the state creation, processing, 
or storage to mimic the experimental procedure. By simu- 
lating smaller or larger imperfections, we can guide the ex- 
perimental focus toward the most critical problems, and es- 
tablish fidelity targets which experimentalists must meet for 
systems to be operationally viable. In network operations, 
the results of a set of experiments run during link initializa- 
tion will be used to track expected behavior of the system. 

Purification takes two imperfect (F < 1.0) Bell pairs 
that span one or more hops and attempts to create one higher- 
fidelity Bell pair with fidelity F\ (F < F' < 1.0). The 
biggest disadvantage of purification is that it is probabilistic, 
requiring two-way classical communication, tying up buffer 
memory resources and reducing throughput. 

2.4 Quantum Repeaters 

The devices that create Bell pairs over a distance are called 
quantum repeaters, building on the concept of teleportation. 
Because teleportation is imperfect, some form of error con- 
trol is necessary. In this paper, we focus on the form of re- 
peater we will call the purify-and-swap approach. Purify- 
and-swap repeaters, as the name suggests, perform two main 
functions: entanglement swapping to lengthen connections 
and purification to manage errors [8. 24 1. Bell pairs are con- 
sumed during the course of teleportation, purification and 
entanglement swapping; thus, the primary job of a repeater 
is to continually produce new ones. From Bell pairs span- 
ning shorter distances, we can create end-to-end Bell pairs, 
which are then given to applications. "Throughput" can be 
defined as Bell pairs per second of a given fidelity, for a 
given path. 

We have defined a protocol stack for purify-and-swap re- 
peaters 11321 . consisting of three layers: Entanglement Con- 
trol (EC) at the link level, and Purification Control (PC) 
and Entanglement Swapping Control (ESC). The latter two 
operate at arbitrary distance. They can be composed recur- 
sively l33l to build high-fidelity Bell pairs in a nested fash- 
ion across any nodes in a network, but doing so requests 
resources which depend critically on the route chosen in the 
network. The operations of these protocol stack elements are 
therefore the ones we vary in the present work. 



As we shall see in Sec. 14.21 the low probability of en- 
tanglement success requires that link-layer operations be ac- 
knowledged. Coupled with very limited buffering resources, 
this results in behavior that differs from most classical net- 
works. Quantum memory is very expensive, simply because 
of the (usually cryogenic) hardware needed to maintain quan- 
tum states for timescales comparable to photon time-of-flight, 
and to connect that hardware to photons. The fidelity of 
data in quantum memory in most technologies decays very 
quickly, adding a hard real-time component to communica- 
tion, or a demand for strong error correction. This decay may 
make congestion control a different, even more severe, prob- 
lem than it is in classical networks. 

For a classical path consisting of a series of identical 
links, throughput (at least in the ideal, sustained case) can 
match the throughput of a link, independent of the number 
of hops. However, the maximum throughput of a symmetric 
chain of purify-and-swap repeaters declines polynomially 
with length, due to the need for additional purification 11241 . 
Thus, it may be desirable to assign a path cost that grows 
more than additively. 

When operated to deliver Bell pairs above a specified 
fidelity threshold (typically related to the fidelity threshold 
for executing quantum error correction), simulations show a 
stair-step decline in throughput versus hop count [32]. The 
total work performed grows linearly as the throughput holds 
steady, then increases suddenly when a "stair" is crossed 
(see Figs. |3]and |4]in SecHJi. Likewise, dealing with non- 
power-of-two numbers of hops will affect throughput and to- 
tal work performed in hard-to-predict fashion, due to changes 
in the swapping and purification patterns. 

These principles have been applied within designs for 
system-area networks, as well as wide-area networks. In- 
terconnects are critical components of quantum computing 
systems [34 35 36 37); even some laboratory-sized experi- 
ments, much like classical multicomputers, will have multi- 
ple levels of interconnect, with different qualities and speeds 
of connection 11381 . These and other large-scale systems will 
use purified Bell states for data movement. Purification forms 
a large chunk of the work in these systems as well as in 
larger networks, and its demands are a commonly-used met- 
ric for evaluating systems. 

Recent results suggest that using quantum error correc- 
tion instead of purification opens up new operational ap- 
proaches more akin to the hop-by-hop behavior of packet- 
switched networks, and may prove to be more efficient in 
many cases 125ll26ll27l . These new approaches reduce the 
demands on memory lifetime and limit the number of re- 
turn messages required. However, the memory resources re- 
quired are substantial, far exceeding what will be practical 
in the near future. The earliest deployments of repeater net- 
works, and perhaps the first commercial products, will likely 
use the purify-and-swap approach. Thus, we reserve analy- 
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sis of our algorithm for other types of networks (and hetero- 
geneous combinations of network types) for future work. 

Although some authors have begun speculating about a 
"quantum Internet", the issues of complex, entangled net- 
works, and heterogeneity at any level above the physical ex- 
change of qubits from one medium to another, are largely 
open problems Q9). 

3 Problems and Solutions 

Our goal for this paper is to develop and analyze a routing 
algorithm for heterogeneous quantum repeater networks, as 
shown in Fig.Q] Our metric for success, therefore, is agree- 
ment between the prospective functioning of a network (rep- 
resented by detailed simulations) and an easy-to-calculate 
algorithmic cost: does our algorithm allow us to make effec- 
tive choices? 

As noted above, to develop a routing algorithm, we need 
a definition for the cost to use a link, a function to calculate 
a path cost based on a set of link costs, and a goal for the 
algorithm (e.g., a metric for deciding if the algorithm meets 
our needs). More precisely, we set as our problems: 

PS.l Choose a goal for the routing algorithm; 

PS.2 for a quantum link, identify the characteristics of in- 
terest for routing, and reduce them to one number or a 
small set of numbers that represent the link cost; and 

PS.3 for a path (an ordered set of links, with associated 
costs), define a function that gives a path cost. 

To solve these problems, we evaluate the following po- 
tential solutions: 

PS.l: Goal As the goal for the routing algorithm itself, 
we choose minimizing work along the path, with attention 
to the secondary goal of selecting the highest-throughput 
single path between the defined communication endpoints, 
measured in Bell pairs per second of the target fidelity. We 
propose a goal for the system of delivering Bell pairs useful 
for teleportation, with a target fidelity of F — 0.98. The phe- 
nomena presented here are independent of the exact value 
chosen, but this value will be adequate for various uses, and 
allows fairly direct comparison to our prior work [t32l . 

Assessing the work for a particular path is not a sim- 
ple problem. Intuitively, we want our measure of work to 
reflect use of some scarce resource. We evaluate two can- 
didates discussed as key functions in repeaters above: total 
measurement operations and total pulse count. 

Without some reasonable idea of global traffic (which 
we don't yet have)Q, evaluation of global success is difficult, 

1 However, for heterogeneous traffic matrices, a gravity model 
might be a reasonable first guess 11391 , as the human patterns driving 
quantum communication are expected to be similar to those for classi- 
cal communication. 



so as a first step we are evaluating the correlation between 
total work on the path, achieved connection throughput and 
our definition of path cost. 

PS.2: Link Cost By analogy with the "transceiver time" 
definition often used with OSPF, we are exploring several 
link cost definitions. Simply assessing the clock speed at 
which laser pulses can be emitted is clearly an inadequate 
link cost metric in a system where the fidelity of the output 
state is important and the probability of success depends on 
the characteristics of the link. Prior to execution of our ex- 
periments, differing intuitions led one author to suggest Bell 
pair generation time and another to suggest the number of 
measurements; eventually we had a list of five candidates: 

Loss the loss in the channel, in decibels; 

InvTrans the inverse of the transmittance of the channel 
(1/T, where the transmittance T is the percentage of 
photons received through the path); 

Pulse the number of laser pulses used to create an entan- 
gled Bell pair of a high fidelity over a single hop, corre- 
sponding to the number of uses of the transmitter (each 
fixed-time); 

Meas the number of measurement operations used to create 
an entangled Bell pair of a high fidelity over a single 
hop (this differs from the above because (a) some pulses 
are discarded rather than accepted when receiving qubit 
resources are busy, and (b) measurements are used in 
entanglement swapping and purification); and 

BellGenT the inverse of the throughput of the single link, 
when run as a single-hop system, measured in seconds 
per Bell pair. 

This list can be divided into two groups: the first two 
candidates are simply physical characteristics of the link that 
can be measured easily, while the latter three require simu- 
lation or monitoring of the link to determine. The first two 
differ by a logarithmic factor; InvTrans corresponds to addi- 
tion while Loss in dB corresponds to multiplication of cost 
when placed in a Dijkstra context. While at first glance it 
may appear desirable to have such an easily-determined link 
cost, by tying the cost so directly to the physical mechanism, 
the definition may not transfer well when links of heteroge- 
neous physical technologies are involved. 

Pulse and Meas seem to correspond most closely to the 
"transceiver time" definition, but BellGenT incorporates sys- 
tem factors and may give more accurate estimates for little 
additional complexity. (Although the technical details are 
very different, these could be considered roughly analogous 
to the raw transceiver rate, the throughput of a flow-controlled 
link, and the throughput of a reliable link-layer protocol, in 
terms of the increasing functionality present.) 

PS.3: Path Cost Function The principle hypothesis of 
our paper is that Dijkstra's algorithm can be used as-is, with 
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an appropriate choice of link cost. More formally, 

C pa th=^C ll ie {P}, (1) 

i 

where {P} is the set of links in a path and c, is the cost for 
link i. When it is necessary to distinguish among the link 
cost candidates, we will refer to them as Dijkstra/BellGenT 
and similarly. 

4 Simulation and Results 

Our goal is to determine the range of conditions under which 
Dijkstra correctly selects the lowest-work path and highest- 
throughput path, and when it selects some reasonable ap- 
proximation. We also wish to articulate the conditions under 
which lowest work and highest throughput are not the same 
path and when the algorithm does not pick the lowest-work 
path. We wish to examine whether we consider those to be 
acceptable cases, or if the algorithm is "failing". This can be 
achieved by creating a set of candidate paths and compar- 
ing the ordering established by Dijkstra/BellGenT and the 
ordering according to simulation of the whole path. 

In this section, we begin by asking a set of questions 
about the behavior of heterogeneous paths (14. 11 1. Next, we 
describe our simulator, proposed hardware configuration and 
single-hop simulation results, which both set the Dijkstra pa- 
rameters and allow us to evaluate our link cost candidates 
( 14.2b . After enumerating a set of interesting path candidates 
( 14.3b . we can answer our questions (14.4b . and use those an- 
swers to solve the research problems posed in Sec. |3~1 ( 14.51 . 

4. 1 The Behavior Questions 

We use the simulator to assess the behavior of systems that 
are too complex to solve analytically and cannot yet be built 
and measured directly. Such simulation results will help to 
guide the development of actual hardware. We can pose a 
series of specific questions that will help us understand how 
heterogeneous paths will behave: 

1 . How does the number of hops in a path affect throughput 
and total work? As a specific case, do the throughputs of 
2" - 1, 2", and 2" + 1-hop paths vary? 

2. How does the number of weak links matter? Does the 
introduction of a single weak link become a bottleneck, 
as in classical networks? Does adding a second or third 
weak link further reduce throughput? 

3. How does the position of weak links in the chain affect 
throughput? Does a weak link at the beginning, in the 
middle, or at the end differ? 

4. Under what circumstances will the path cost mis-order 
candidate paths with respect to throughput or total work? 



For classical systems, we know the answers: although 
cost increases, (1) the length of the path (in hops) has no 
effect on the (theoretical) throughput, and (2) the (theoreti- 
cal) throughput of a path is capped by the throughput of the 
bottleneck link(s). Question (3) is answered in the negative: 
ordering does not (or should not) matter, and ordering of 
path segments should always agree with ordering of the full 
paths. Question (4) is open-ended, but there are known cases 
where a lower-performance path can be selected in the ser- 
vice of a larger global goal. The answers to these questions 
are presented in Sec. 14.51 



4.2 Simulated Hardware and Link Costs 

All of the simulations presented in this paper model the 
aubus physical entanglement mechanism 14011411 . Qubus uses 
a strong light pulse to entangle two distant, static qubit mem- 
ories through a fiber or other waveguide. Alternative ap- 
proaches use single photons l42l or a small number of pho- 
tons 0431 . In the qubus system, a strong laser pulse first in- 
teracts with a physical qubit at the transmitter, such as an 
atom or a quantum dot held in a cavity, resulting in a small, 
nonlinear shift in the state of the laser pulse, depending on 
the state of the qubit. The pulse is sent through a waveguide 
to the receiver, where it undergoes a similar interaction with 
another qubit. The pulse is then measured using a technique 
known as homodyne measurement, giving a classical result. 
About 40% of the time, we get a result that tells us only the 
parity of the two qubits and leaves them entangled. 

We use the same simulator as in our prior work, with ex- 
tensions to support the heterogeneous paths [32 40 43 1. The 
simulator was developed for modeling the quantum-level 
behavior of a cavity QED system for qubus and two other 
physical layer candidates, and uses the purification mecha- 
nism of Dehaene et al. fi4l . The necessary classical mes- 
saging is carefully modeled. The simulator consists of about 
11,000 commented lines of C++, and the production runs 
of the simulations presented here consumed about one hun- 
dred hours of CPU time on 2.2GHz AMD Opteron CPUs. 
The simulator is based on well-understood physical equa- 
tions and is experimentally validated at the lowest levels; 
as experience with larger-scale quantum networks develops, 
we plan to continue tracking the agreement of theory, sim- 
ulation, and experiment. We believe the prospective agree- 
ment between a real-world network and our simulations will 
be more than adequate for the purposes of this paper, but the 
detailed arguments are largely about physics, and are beyond 
the scope of this paper. 

Figure [3] shows the results of simulating a single qubus 
hop, with parameters as in Table \T\ additional details of the 
hardware configuration are the same as in Refs. 132 4011431. 
The original qubus mechanism is very sensitive to loss, but 
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works well in low-loss situations. This form of qubus re- 
peater fails to work for losses greater than about 5.5dB from 
transmitter qubit to receiver qubit, limiting hop length to 
about 30km over high-quality optical fiber at telecom wave- 
lengths. Other types of physical link, including variants of 
the qubus mechanism B31 . will work over longer distances 
and with higher fidelity. By choosing to simulate basic qubus 
links, we can see very clearly the impact of purification and 
low-fidelity entanglement. As future work, we plan to con- 
firm the behavior of Dijkstra with other physical link types. 

For the multi-hop simulations, we choose four specific 
points as example links, as shown in Table [2] The four link 
types chosen are marked in the figure with the corresponding 
symbols. The terms "standard", "good", "fair", and "poor" 
are relative to this simulation only, not indicative of all pos- 
sible physical quantum link types. As can be seen from this 
table and Figure[3] even small differences in loss have a large 
and uneven impact on throughput, suggesting that not only 
is the utility of the proposed link costs Loss and Inv Trans 
rather technology-specific, but even in the isolated case of 
qubus they may be poor metrics. 




3.4 3.6 3.8 4.0 4.2 4.4 

Loss (dB) 



Fig. 3 Single-hop simulated throughput versus loss using the parame- 
ters in Table \l] used to define linkcost. The four link types chosen for 
simulation of more complex paths are marked on the black/blue curve 
with the corresponding symbols. The stair-step behavior is due to the 
increasing number of purification rounds as the initial fidelity declines 
due to increasing loss. The details of this curve are very specific to the 
qubus technology we simulate, but the principles are general. 



Although Bell pairs are symmetric and can be used to 
teleport data in either direction, at the physical level, links 
are naturally dependent on a transmitter and a receiver. Re- 
ceivers dynamically select a qubit to attempt to entangle, 



Table 1 Link cost simulation parameters 



Qubus Quantum Repeater 


Number of qubits 


25 transmitter, 


per repeater link connection 


25 receiver 


Number of qubits 
teleported (length of simulation) 
Final target fidelity 


200 qubits 
F > 0.98 


Optical Fiber 


Length 


20 km 


Signal loss 


3.4-4.4 dB/20 km 



based on availability at the time of pulse reception, much 
as any network interface chooses a buffer in which to place 
an incoming packet. The inherent round-trip latency in re- 
ceiving acknowledgement of entanglement success or fail- 
ure (the Entanglement Control (EC) protocol layer) means 
that the qubits at the transmitter must spend long periods of 
time in which they are nominally marked as busy, but no 
useful forward progress can be be made. Thus, links with 
the same physical loss but different numbers of qubits at the 
transmitter and receiver will behave differently, as shown in 
Table[3] 

Although the throughput varies by a factor of two for 
different configurations, the number of pulses and number 
of measurements to teleport 200 qubits does not vary sig- 
nificantly except for the clearly misconfigured 100/25 case. 
For the conditions simulated, about 450 pulses and 700 mea- 
surements (including basic entanglement, swapping and pu- 
rification) are required to teleport one qubit. Although these 
two measures are arguably more direct representations of 
cost, this important throughput difference affects our notion 
of a preferred path. 

Reasoning about the behavior of only a single hop and 
confirming the results via simulation allow us easily to elim- 
inate two of our prospective link cost candidates (Loss and 
InvTrans). Two others (Pulse and Meas) remain good mea- 
sures of total work, but the results shown in Tab. [3]cast doubt 
on their viability as candidates for link cost when the goal is 
to achieve high throughput. Thus, we settle on BellGenT as 
our link cost metric. 

4.3 Simulated Path Candidates 

The candidate paths we choose to examine are designed to 
answer the questions in Sec. 14.11 We simulated forty-six 
paths of one to nine hops in various patterns, as well as all 
4 4 = 256 four-hop combinations for the four chosen link 
types. The forty-six paths are enumerated along the bottom 
of Fig. [4] We have simulated many more paths with a va- 
riety of link conditions over the course of this experiment, 
notably very long, homogeneous paths (up to 2,048 hops). 
The findings from other simulations do not contradict the 
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Table 2 The link configurations and five candidates for link cost. The cost candidates to the left of the first double line are physical characteristics 
of the channel, while the others are the results of single-hop simulations. Simulation parameters are as in Table[T] 



Link 


Loss 


InvTrans 


Pulse 


per-t'port 


(norm) 


Meas 


per-t'port 


(norm) 


T'put 


BellGenT 


Standard (□) 


3.4dB 


2.19 


90441 


452 


1 


140519 


702 


1 


217.7 


1 


Good (O) 


3.7dB 


2.34 


163628 


818 


1.80 


254691 


1237 


1.76 


118.4 


1.83 


Fair (A) 


4.0dB 


2.51 


258852 


1294 


2.86 


404117 


2020 


2.87 


74.3 


2.93 


Poor ( x ) 


4.3dB 


2.69 


606278 


3031 


6.70 


945247 


4276 


6.72 


33.1 


6.57 



Table 3 Throughput depends on the number of available qubits at both 
transmitter and receiver. Due to the round-trip latency in the Entangle- 
ment Control (EC) protocol for acknowledging successful entangle- 
ment, having more qubits at the transmitter boosts throughput. Simu- 
lations are for a single hop with parameters as in Table [T] Throughput 
and confidence interval (std. dev.) are determined by a linear fit to tele- 
portation completion times. 



Xmtr 


Rcvr 


Throughput 


Pulses 


Meas 


25 


25 


237 ± 8 


80587 


125905 


25 


50 


213 ± 8 


92071 


143582 


50 


25 


436 ± 14 


99708 


142723 


50 


50 


456 ± 14 


89533 


139506 


100 


25 


462 ± 22 


199240 


138543 


100 


50 


956 ± 31 


92221 


137828 


100 


100 


984 ± 32 


88694 


138145 



results presented in this paper, which were chosen to clearly 
show the effects of interest. 



4.4 Answering Our Behavior Questions 

Figures|4][5]and[6]show our simulation results for the paths. 
In Fig. 2J the throughput for each specific path can be seen, 
as well as the two measures of total work, Pulse and Meas. In 
Fig. |5j throughput is plotted against the Dijkstra-calculated 
path cost using BellGenT as our cost metric, and in Fig. [6] 
the total work measures are plotted against calculated path 
cost. Figures Q and [8] plot the results for all 256 four-hop 
paths we simulated. 

Perhaps the most important factor to note about the be- 
havior we see is the discrete nature of changes to throughput, 
due to the discrete nature of purification and entanglement 
swapping and our choice to establish a particular threshold 
for final acceptance of an end-to-end Bell pair. In many cir- 
cumstances, minor changes force an additional round of pu- 
rification, generally causing a 50% reduction in throughput. 

As the number of hops (Q.[U) increases, our results show 
a stair-step phenomenon in throughput, but not necessarily 
at powers of two: the decline of fidelity resulting in addi- 
tional rounds of purification does not move in concert with 
entanglement swapping. For Standard □ paths, we see steps 
at the 2nd and 7th hops. 

Examining the symbols in Fig. [5] shows clearly the ex- 
istence of a bottleneck link phenomenon (Q. [2]): with one 
low-quality link in the path, the quality of the other links is 
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Fig. 4 Total number of pulses (A) and measurements (+) for forty- 
six of our candidate paths (right scale). The paths vary in length from 
one to nine hops. They are ordered left to right according to ascending 
throughput, plotted using bars (left scale). The legend below the graph 
shows the individual path configurations; □, 0> A and x represent 
our standard, good, fair, and poor links, respectively. The stair-step be- 
havior reflects increasing numbers of rounds of purification. 



almost irrelevant, as can be seen by the clustering of each 
type of data point (e.g., A at about 30 Bell pairs/second in 
Fig. [5] and the concentration of x links at the left edge of 
Fig. 0}. The various paths can largely, though not entirely, 
be grouped according to the throughput of the slowest link 
in the path. The most eye-catching anomalies in Fig. [5] are 
the single 0> A and x marks above and to the right of their 
respective clusters. These are the single-hop paths of the cor- 
responding link type, indicating that we do not have a pure 
bottleneck phenomenon. The bottleneck plus the polynomial 
decline in performance as the number of hops grows work 
in combination determine to the path performance. In most 
cases for the longer paths, adding a second link of the same 
quality as the bottleneck link does not result in a statistically 
significant reduction in throughput, but optimization of the 
path becomes more difficult. In a few cases, our optimization 
fails to find an acceptable pattern, and an additional round of 
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Fig. 5 Throughput versus BellGenT path cost for forty-six of our can- 
didate paths. Each path is represented by the symbol for the weakest 
type of link in the path. The clustering of each type of data point shows 
clearly that throughput is limited by the bottleneck link. The length of 
the vertical bar (mostly contained within each symbol) shows the std. 
dev. of the throughput. 



purification becomes necessary. In particular, four-hop paths 
with O as m e bottleneck link(s) are split into two plateaus at 
around 80 Bell pairs/second and 50 Bell pairs/second. How- 
ever, the total work shows strong correlation even for these 
cases. 

Despite this general bottleneck behavior, the weak link 
position (Q.|3]l is a more subtle one. Among our simulations, 
we found a single case where the bottleneck position pro- 
duced different results. With one Good link and three Stan- 
dard ones, having the Good link at the left end of the path re- 
sulted in 53 Bell pairs/second, whereas the other three paths 
produced 82 to 87 Bell pairs/second. The low-throughput 
case required an extra round of purification before one en- 
tanglement swapping operation. As the path hovers near a 
threshold demanding an extra round of purification, opti- 
mization of the path usage becomes more difficult; three 
paths successfully did so, while the fourth didn't. Thus, we 
must answer that the weak link position may have an impact 
on throughput, depending on our ability to effectively use 
the path. 

Dijkstra/BellGenT occasionally mis-orders (Q. [4]i pairs 
of path candidates with respect to throughput. Comparing 
the 256 four-hop paths we simulated, there are 32,640 possi- 
ble pairs, of which 1,230 had the same path cost. Of the pairs 
with different costs, 82.6% of the time Dijkstra/BellGenT 
chooses the higher-throughput path of the pair (the "cor- 
rect" choice), and 17.4% of the time it chooses the lower- 
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Fig. 6 Total work, measured in Pulses (A) and Measurements (+), 
versus BellGenT path cost for forty-six of our candidate paths. The 
coefficient of determination of each linear fit is 0.88, showing that our 
path cost is a strong predictor of total work. 
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Fig. 7 Throughput versus BellGenT path cost for all 256 four-hop can- 
didate paths. Each path is represented by the symbol for the weakest 
type of link in the path. The clustering of each type of data point shows 
clearly that throughput is limited by the bottleneck link. 



throughput path (the "incorrect" choice). In only 5% of those 
mis-ordering cases was the difference in throughput more 
than 10%. For the 46 variable-length paths, the rate of "cor- 
rect" choices was similar, 81.6%, but the impact of those 
choices is higher, with a throughput penalty of 25% or more 
in almost half of the "incorrect" cases. 
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Fig. 8 Total work, measured in pulses (A) and measurements (+), 
versus BellGenT path cost for all 256 four-hop candidate paths, with 
linear fits. The coefficient of determination for the number of pulses is 
0.81, and for the number of measurements is 0.99. 



The most remarkable case of incorrect ordering we found 
is the pair of paths □ x □□□□□□ (throughput: 16.0 ±0.51 
cost: 13.57) and AAAADDDD (throughput: 30.7 ± 0.77 
cost: 15.72). The higher-cost path has twice the throughput 
of the lower-cost one. Examining Fig.|4]shows that the total 
work (in pulses or measurements) is quite similar. This iso- 
lated case would suggest that the bottleneck Poor link might 
warrant even a higher link cost than 6.57, however, in all of 
the other cases we examined, it has worked well. 

The performance and work are both dominated by pu- 
rification costs. The use of a given path must be optimized 
as a whole; our current solution is brute force, trying a large 
number of possibilities. In a large majority of cases, opti- 
mization of that process for a path is straightforward and 
robust. However, a noticeable minority of cases require deli- 
cate adjustments to the entanglement swapping settings. Fur- 
ther automated optimization of this process may result in 
both better performance and agreement between path cost 
and work. 

For total work, we find a coefficient of determination 
of 0.88 for both pulses and measurements for the forty-six 
variable-length paths. For the 256 four-hop paths, we find a 
coefficient of 0.8 1 for pulses and 0.99 for measurements. 

We can summarize the behavior as follows: Performance 
is determined by the number of rounds of purification used 
anywhere on the path; the number of rounds is dominated 
by the bottleneck link. Work, however, is spread across the 
entire path in rough proportion to relative link quality. 



4.5 Solving Our Problems 

With the simulator results and answers to our questions in 
hand, we are now prepared to assess the solutions we pro- 
posed in Sec. [3] and determine whether Dijkstra/BellGenT 
meets our goal of developing an acceptable routing algo- 
rithm for quantum repeater networks. 

We have seen that the highest-throughput path is also 
strongly correlated to measures of work, including total pulse 
count and number of measurements. In some pathological 
cases, number of measurements is a better metric, and across 
a broad range of cases it matches well with our chosen link 
cost below. Our solution to Path Selection Problem PS.l 
is therefore: to minimize total work along the path, using 
number of measurements as the metric. 

We saw in Sec. 14.21 that defining BellGenT as our link 
cost metric is likely to better suit our purposes than Pulse, 
Meas, Loss, or InvTrans, allowing us to propose a solution 
to Problem PS.2: seconds per Bell pair. 

Our simulation results presented in the previous sub- 
sections indicate clearly that we can usually, if not always, 
correctly predict the highest-throughput, lowest-work path 
from among a set of candidates. Dijkstra/BellGenT solves 
Problem PS.3 well: path cost is the scalar sum of link costs. 

Intuitively, we can see why this works well, as confirmed 
by the simulations: BellGenT is directly related to the num- 
ber of purification rounds needed, and hence the amount 
of work on a link. Despite the concerns we expressed ear- 
lier about the nonlinear amount of work as the number of 
hops grows, the predictive ability for ordering paths remains 
strong. Thus, we can assert that even as the use of complex 
repeater networks evolves, and various traffic patterns arise, 
Dijkstra/BellGenT likely will remain an effective choice. 



5 Conclusion 

Real-world deployment of networks of quantum repeaters 
will inevitably be physically heterogeneous, with complex 
topologies and high- and low-quality links and many pos- 
sible paths through the networks, rather than an idealized, 
homogeneous power-of-two number-of-hops. 

In this paper, we have investigated critical problems in 
the use of purify-and-swap repeater networks. Our focus 
has been on path selection, and the need for a routing al- 
gorithm. An acceptable routing algorithm must be an easy- 
to-calculate metric that reliably chooses a reasonable, if not 
optimal, path. 

Our results show that, despite many important differ- 
ences, quantum repeater networks behave similarly to classi- 
cal networks in useful ways, but the classical principles can- 
not be applied without thought. Via detailed physical simu- 
lation of both the physics and the classical messaging proto- 
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cols, we have investigated several variants and explored the 
range of conditions under which these principles apply. 

We can predict the throughput of a connection based pri- 
marily on the bottleneck link in the path, while the total 
work, in number of operations performed (pulses or mea- 
surements), increases with the addition of other, non-bottleneck 
links, much as in a classical network path. Applying a form 
of Dijkstra's algorithm with the inverse throughput of each 
hop as the link cost (Dijkstra/BellGenT) results in strong 
correlation between our easily-calculated path cost and ac- 
tual throughput, and between cost and total work. This is 
achieved with reasonable computational effort, allowing us 
to recommend the use of Dijkstra. 

In future work, we expect to confirm the general be- 
havior of Dijkstra with other physical-layer repeater types, 
adding single-photon P31 and low-photon number B31 to 
the qubus systems explored here. Our simulator is capable 
of modeling the important factor of finite quantum memory 
lifetimes ||46| . but the sheer additional combinatoric com- 
plexity that would have come from including variation of 
this parameter in both the simulations and algorithmic argu- 
ments in this paper prevents us from presenting these results. 

Perhaps the most important open question is whether 
these results apply to error-corrected, rather than purified, 
repeater networks H25 112611271 . Demonstrating the applicabil- 
ity of Dijkstra to both purify-and-swap and error-corrected 
repeater networks would be a strong indicator of the univer- 
sality of our results. We expect that initial network demon- 
strations will be purify-and-swap, moving toward error-corrected 
but perhaps with a long period of coexistence. 

The results presented here, concerning routing within a 
network of uniform technology but disparate operating con- 
ditions, are part of a larger program attempting to unify var- 
ious quantum networking approaches within a single frame- 
work, covering inter-networking between heterogeneous phys- 
ical technologies as well as radically different error correc- 
tion approaches, routing, multiplexing approaches and qual- 
ity of service, and distributed applications. With these re- 
sults, we expect to take advantage of the knowledge accu- 
mulated in half a century of classical networking research 
and operation to create a strong quantum network architec- 
ture, guide the experimental focus for the development of 
repeaters, implement real-world quantum protocols and ac- 
celerate deployment of quantum networks. 
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